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摘要 :【 目的] 拟 步 甲 科 昆虫 小 胸 警 甲 Microdera punctipennis 是 分 布 于 中 国 新 疆 古 尔 班 通 古 特 沙漠 
的 特有 物种 ,具有 很 强 的 耐寒 性 ,但 其 低温 响应 的 分 子 机 制 尚 不 明确 。 本 研究 旨 在 利用 转录 组 测序 
技术 丰富 小 胸 束 甲 的 已 知 基 因 信 息 ,分析 在 短 时 低温 胁迫 下 ,上 调 表达 基因 的 主要 类 和 群 及 参与 的 主 
要 代谢 通路 。【 方 法】 利用 Trinity 软件 对 分 别 在 4% 低温 和 25% (对 照 ) 下 处 理 3h 的 小 胸 整 甲 成 
X RNA-Seq 数据 进行 从 头 组 装 ,利用 Blas2go 软件 进行 基因 注释 。 通 过 DESeq 和 GFOLD 软件 分 
析 差 异 表达 基因 ,并 对 上 调 表达 基因 进行 GO(Gene Ontology) fe KEGG ( Kyoto Encyclopedia of Genes 
and Genomes) 代谢 途径 富 集 分 析 。【 结果】 测序 过 滤 后 得 到 68 165 001 个 有 效 读 序 ,51 712 个 平均 
长 度 为 984 bp 的 非 宛 余 基 因 , 其 中 29 925 个 基因 ( 约 57.87% ) 有 同 源 序 列 (E-value «10 7) ,与 拟 
步 甲 科 的 赤 拟 谷 盗 Tribolium castaneum 相应 基因 相似 性 最 高 ( 核 苷 酸 序列 一 致 性 为 72.75% ) 。 低 
温 响 应 上 调 表达 基因 有 514 个 , 富 集 到 35 个 GO 类 群 ,18 个 KEGG 途径 。 其 中 胁迫 响应 类 群生 物 
调节 类 群 和 免疫 系统 过 程 类 群 都 占有 重要 比例 , 嗓 叭 代谢 、 硫 胺 素 代谢 、 糖 酵 解 / 糖 异 生 等 代谢 通路 
显著 富 集 。 在 胁 连 响应 类 群 中 , 热 激 蛋白 Hsp90 基因 、 超 氧化 物 歧化 酶 (SOD ) 基 因 和 钙 联 接 蛋 和 白 
Calnexin 基因 等 8 个 非 生 物 胁迫 相关 基因 显著 上 调 表达 。【 结论 ] 荒漠 昆虫 小 胸 鳌 甲 低温 转录 组 揭 
示 , 在 4% 冷 驯化 过 程 中 ,代谢 过 程 、 胁 迫 响 应 生物 调节 和 免疫 系统 等 生物 学 过 程 相关 基因 表达 显 
著 上 调 , 其 中 几 个 非 生物 胁迫 响应 基因 提示 小 胸 鳌 甲 可 能 从 分 子 伴侣 、 细 胞 周期 阻 滞 、 线 粒 体 稳定 
性 、 抗 氧化 胁迫 等 多 方面 应 对 低温 胁迫 。KRGG 富 集 分 析 所 揭示 的 小 胸 鳌 甲 在 低温 下 的 转录 组 代 
谢 通 路 与 其 他 物种 有 较 大 差异 。 研 究 结果 为 后 续 生物 信息 学 分 析 及 小 胸 鳌 甲 耐 寒 关 键 基 因 的 发 据 
及 耐寒 机 制 的 揭示 提供 了 基础 数据 。 
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Transcriptomic analysis of the desert beetle Microdera punctipennis 


( Coleoptera: Tenebrionidae) in response to short-term cold stress 
Kuerban TUSONG, LU Xue-Ying, LIU Xiao-Ning, MA Ji" ( Xinjiang Key Laboratory of Biological 
Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, 
Urumqi 830046 , China) 

Abstract: [ Aim] Tenebrionid insect Microdera punctipennis is an endemic species in Gurbantünggüt 
Desert, Xinjiang, China. It has very strong cold hardiness, but the molecular mechanisms by which M. 
punctipennis responses to low temperature are still unclear. This study aims to enrich the genetic 
information of M. punctipennis , and to analyze the functional groups of the most up-regulated genes under 


short-term cold stress and the pathways which they may involve in. [ Methods] Transcriptomic sequencing 
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reads from M. punctipennis adults exposed to 4°C as low temperature treatment and 25°C as the control, 
respectively, for 3 h were de novo assembled by using Trinity software. Gene annotation was conducted by 
using Blast2 go software. The differentially expressed genes were screened by using DESeq and GFOLD 
softwares. For the up-regulated genes, GO term enrichment and KEGG metabolic pathway analysis were 
performed using Blast2go software. [Results] In total, 51 712 non-redundant transcripts with the average 
length of 984 bp were assembled, of which the expression levels of 514 genes were up-regulated and 91 
genes down-regulated in response to cold stress. The up-regulated ones were enriched to 35 GO terms and 
18 KEGG pathways, showing that stress-related GO terms, biological regulation and immune system are 
enriched significantly. Pathways such as purine metabolism, thiamine metabolism, and glycolysis/ 
gluconeogenesis are significantly enriched. In the stress responsive GO term, the expression levels of 
eight abiotic related genes, including heat shock protein ( Hsp90) gene, superoxide dismutase ( SOD) 
gene, calnexin gene, etc., were up-regulated significantly. [Conclusion] The transcriptomic data of the 
desert beetle M. punctipennis reveal that the expressions of genes related to the biological processes like 
metabolic process, stress response, biological regulation and immune system are significantly up-regulated 
during cold acclimation at 4°C. The metabolic pathways revealed by KEGG enrichment show that low 
temperature transcriptome of M. punctipennis is much different from those of other species. Our results 


provide basic data for successive bioinformaties analysis and mining the key genes in cold tolerance for 


M. punctipennis as well as for elucidation of the related molecular mechanisms. 
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昆虫 占 地 球 动物 种 类 的 75% 以 上 ,存在 于 几乎 
所 有 的 生态 系统 中 , 并 发 挥 着 重要 的 生态 作用 
(Chown and Nicolson, 2004; Harrison et al., 2012; 
MacMillan et al., 2015) 。 温 度 通过 约束 昆虫 的 地 理 
分 布 和 季节 性 活动 ( Harrison et al., 2012) , 直接 或 
间接 地 影响 昆虫 的 扩散 和 侵害 性 、 天 敌 昆 虫 引 入 的 
成 功 性 , 及 原生 昆虫 种 群 的 动态 性 (Bale and 
Hayward, 2010)。 极 地 和 高 山地 区 的 大 部 分 昆虫 能 
够 适应 低温 及 冬季 的 环境 胁迫 而 生存 (Leather et 
al., 1993) ,荒漠 昆虫 也 经 常 遭 遇 潜 在 的 致命 寒冷 气 
Ik (Sømme, 1995) 。 因 此 ,低温 生物 学 是 昆虫 适应 
性 的 重要 组 成 部 分 ,也 是 昆虫 分 布 的 重要 决定 因素 
之 一 (Andersen et al., 2015) 。 

在 低温 条 件 下 ,昆虫 会 发 生 一 系列 生理 生化 改 
变 (Denlinger and Lee, 2010) 。 近 年 来 ,从 基因 转录 
组 水 平 研究 昆虫 对 低温 的 适应 机 制 进展 迅速 。 例 
如 ,利用 抑制 消减 杂交 技术 (SSH ) 研究 尖 音 库 蚊 
Culex pipiens 冬季 清 育 成 虫 (Robich et al., 2007) 分 
析 冷 驯化 及 冷 休 克 复 苏 后 的 黑 腹 果 晶 Drosophila 
melanogaster 转录 组 (Denlinger and Lee, 2010) 以 及 
红 尾 肉 量 Sarcophaga crassipalpis 转录 组 (Rinehart et 
al., 2007) ,发现 了 几 个 与 冷 应 激 相关 的 基因 上 调 表 
达 。 利 用 基因 芯片 技术 分 析 0% 处 理 2 h 并 复苏 30 
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min 的 果 蝇 ,发 现 膜 蛋白 胁迫 响应 蛋白 和 线粒体 蛋 
白 均 有 差异 表达 ; 7 0. 5C 处 理 不 同时 间 的 果 蝇 转录 
组 ,基因 的 表达 有 很 大 差异 (Qin et al., 2005; Zhang 
et al., 2011), TRR! Sarcophaga bullata 的 低温 
转录 组 ,发 现 快速 冷 驯化 对 其 转录 组 表达 谱 的 影响 
不 大 (Teets et al., 2012) 。 

RNA-Seq 技术 的 广泛 运用 , 极 大 地 促进 了 非 模 
式 生物 的 功能 基因 组 分 析 , 例 如 对 竹 节 虫 Micrarchus 
sp. ,/|NZ£ Tomicus yunnanensis 7H Sf HR, Cryptolaemus 
montrouzieri ( Zhu et al., 2012; Zhang et al., 2012; 
Dunning et al., 2013) m c SRR XI H R 
Ericerus pela 的 过 冷却 点 转录 组 研究 发 现 有 1 124 条 
基因 上 调 表达 ,1 262 条 基因 下 调 表达 ,这 些 差异 性 
表达 基因 分 别 富 集 到 蔗糖 代谢 . 半 乳 糖 代 谢 、 多 种 氮 
基 酸 代谢 等 途径 (Yu et al., 2016) 。 

7f] FH. Microdera punctipennis ( 8$ 38 H : $127 
甲 科 ) 是 分 布 于 新 疆 古 尔 班 通 古 特 沙漠 的 特有 物种 
CLA $E SE, 2005) 。 该 地 区 冬季 最 低 气 温 可 降 至 
-41. 5*C (Zhou et al., 2009) ,. /J ffs] Ex FH AC SE TA 
夜 出 来 躲避 干旱 和 高 温 , 并 大 量 合成 抗 冻 蛋 白 过 冬 ， 
其 成 虫 的 过 冷却 点 在 冬季 可 低 至 -19.6%C (Hou et 
al., 2010) ,其 耐 受 荒漠 极端 环境 的 潜在 分 子 机 制 还 
需 从 基因 转录 组 水 平 进行 深入 研究 。 之 前 我 们 利用 
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RNA-Seq 技术 Illumina Solexa 平台 对 4% 低温 短 时 
胁迫 及 25C 对 照 组 的 小 胸 整 甲 成 虫 进行 了 转录 组 
测序 ,由 于 拼接 组 装 获得 144 777 条 基因 (Lu et al., 
2014b) ,显著 多 于 已 报道 的 昆虫 转录 组 数据 ,可 能 
存在 大 量 元 余 转 录 本 , 故 本 研究 对 测序 数据 进行 了 
重新 组 装 , 共 获 得 51 712 条 基因 ,有 效 地 剔除 了 元 
余 转 录 本 ,并 对 4Y 低温 上 调 差 异 表达 基因 进行 了 
功能 分 类 ,以 期 发 掘 与 小 胸 获 甲 低温 适应 机 制 相关 
的 关键 基因 和 低温 相关 的 代谢 通路 ,为 从 分 子 水 平 
开展 小 胸 整 甲 低温 适应 机 制 研 究竟 定 基础 。 
































1 材料 与 方法 


1.1 试 虫 来 源 及 处 理 

^g E FP p Ha F 2011 年 10 月 2 日 采 自 新 疆 古 
尔 班 通 古 特 沙漠 边缘 地 带 (44°24'N，87°51'F, 海 拔 
444 m) ,在 室内 25 人 饲养 一 周 后 选取 个 体 大 小 相近 
的 10 头 成 虫 用 于 实验 ,随机 将 5 ARR BET ACE 
为 处 理 组 ,剩余 5 KEF 25 C 作为 对 照 组 ,处 理 3 h 
后 迅速 将 处 理 组 和 对 照 组 试 虫 置 于 液 所 中 , 再 转移 
至 -80% 冰箱 以 备 RNA 提取 。 
1.2 cDNA 文库 构建 和 Tllumina 测序 

从 两 组 冻 存 的 小 胸 警 甲 试 虫 中 各 随机 取出 3 
头 ,用 TRIzol 试剂 提取 总 RNA。 总 RNA 的 提取 、 质 
量 分 析 、cDNA 文库 的 建立 及 测序 工作 由 生 工 生物 工 
程 (上 海 ) 有 限 公司 完成 。 采 用 Humina HiSeq™ 2000 
(Illumina, San Diego, CA, USA) 系统 进行 测序 。 
1.3 序列 拼接 及 元 余 转录 本 去 除 

利用 Nesoni tools 程序 对 原始 测序 数据 进行 质 
量 预 处 理 。 质 量 闵 值 为 20( 最 大 错误 率 1% ) ,长 度 
阅 值 为 24 bp ,切除 读 序 中 含 N 的 部 分 。 使 用 Trinity 
v2.0.6 对 得 到 的 有 效 读 序 (clean reads ) 进行 paired- 
end 拼接 。 由 于 利用 Trinity 组 装 过 程 中 会 存在 单 核 
苷 酸 多 态 性 和 一 些小 插入 或 者 缺失 所 造成 的 元 余 
( Wang et al., 2013) , 故 使 用 CD-HIT-EST 软件 删除 
完全 被 长 转录 本 履 盖 且 一 致 性 (identity ) 高 于 9996 
的 短 宛 余 转 录 本 (基因 ) ,得 到 非 元 余 基 因 。 
1.4 生物 信息 学 分 析 

使 用 Blast2go 软件 将 拼接 得 到 的 非 元 余 基因 与 
NCBI 的 NR 和 蛋白 质数 据 库 (最 后 升级 时 间 2015 年 5 
月 ) 进行 BLASTX 注释 (E-value < 10 7) ; Ñ go 
mapping, go annotation, interpro scan 和 merging of 
interpro GOs to annotation 等 步骤 寻找 与 某 些 结构 域 
寺 征 序列 同 源 的 基因 进行 CO (gene ontology) 注 

































































释 ,并 进行 GO term 富 集 分 析 。 利 用 GO 数据 库 提 
fit fj GO-Enzyme Code 文件 和 KEGG ( Kyoto 
Encyclopedia of Genes and Genomes ) 数据 库 提 供 的 
代谢 通路 信息 ,对 参与 代谢 通路 的 序列 的 富 集 程度 
进行 分 类 。 
1.5 小 胸 警 甲 低温 转录 组 差异 表达 分 析 

利用 Bowtie2 软件 对 有 效 读 序 和 非 元 余 基 因 进 
行 mapping, 得 到 alignment 文件 (Sam file) ,通过 
samtools 处 理 得 到 各 转录 本 在 不 同 处 理 条 件 下 的 读 
数 (readcount) 。 把 readcount 文件 作为 收入 文件 ,在 
R 环境 使 用 bioconductor 的 DESeq 工具 分 析 低 温 处 
理 数 据 的 差异 表达 转录 本 , 将 错误 发 现 率 (false 
discovery rate, FDR) 分 别 设 为 小 于 0.05 和 0.1, 获 
得 差异 表达 基因 的 倍数 均 三 3.6 以 此 为 基础 ,将 差 
异 表达 倍数 4 VON BJ RI GFOLD 软件 (Feng et al., 
2012 ) 对 小 胸 警 甲 转录 组 数据 再 次 进行 差异 表达 分 析 。 
取 两 种 方法 的 并 集 作为 第 选 出 的 差异 表达 基因 。 

















2 结果 


2.1 低 水 平 读 序 的 去 除 及 从 头 组 装 结果 

对 4% 处 理 组 和 25%C 对 照 组 小 胸 鳖 甲 成 虫 
cDNA 文库 采用 Ilumina HiSeq” 2000 测序 , 共 获 得 
84 944 120 条 原始 读 序 ,经 过 质量 预 处 理 后 最 终 获 
得 平均 长 度 为 80.6 bp 的 有 效 读 序 68 165 001 条 。 
利用 Trinity v2. 0. 6 对 有 效 读 序 进行 组 装 , 初 步 得 到 
56 130 条 基因 ,去 除 见 余 序列 后 最 终 得 到 51 712 条 
非 元 余 基因 ( 表 1) ,可 对 转录 组 质量 进行 评价 。 
2.2 功能 注释 及 序列 特征 

将 全 部 51 712 个 基因 与 NR 数据 库 进 行 序列 同 
源 性 比 对 ,共有 29 925 个 ( 占 57.87% ) 基因 与 其 他 
生物 的 已 知 蛋 白 具 有 不 同 程度 的 同 源 性 。 在 已 注释 
序列 中 ,完全 匹配 的 有 7 265 个 (E-value 20) , 占 注 
释 序列 的 24.28% ; t FE VU ig ( E-value « 1. OE — 50) 
的 有 9 495 A ( 31. 7396) ; 其 余 13 161 个 ( 占 
43.9896) E fiy 1.0E-50 -1.0E-5 (图 1: A) 

从 所 注释 到 的 物种 分 布 来 看 ,有 21 769 个 ( 占 
72.75% ) 基因 与 拟 步 甲 科 模式 昆虫 赤 拟 谷 盗 T. 
castaneum 同 源 ;2 182 个 ( 占 7.29% ) 基 因 与 节肢 动 
物 寄 生 虫 Gregarina niphandrodes 同 源 ;360 个 ( 占 
1.296) 基因 与 中 欧 山 松 大 小 喜 Dendroctonus 
ponderosae 同 源 ;237 个 ( 占 0.79% ) 基因 与 黄粉 虫 
Tenebrio molitor 同 源 , 与 其 他 物种 同 源 的 基因 均 少 
于 200 个 (图 1: B), 
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Table1 Statistics of the sequencing results and de novo assembly of Microdera punctipennis tanscriptome 
低温 处 理 组 对 照 组 合计 
Cold treatment group Control group Total 
读 序 Reads 
总 原始 读 序 Total raw reads 36 786 116 48 158 004 84 944 120 
原始 读 序 平均 长 度 Average input length ( nt) 100 100 
有 效 读 序 平均 长 度 Average output length (nt) 79.53 81.7 
总 有 效 读 序 Total clean reads 29 066 064 39 098 937 68 165 001 
有 效 读 序 Q20 百分比 Q20 percentage 79. 0196 81.1996 
组 装 基因 Assembled unigenes 
总 数 Total number 39 712 42 409 56 130 
非 兄 余 基因 Filtered unigenes 
总 数 Total number 37 086 40 032 5] 712 
最 长 基因 Maximum length ( nt) 27 511 12 194 30 733 
最 短 基 因 Minimum length ( nt) 224 224 224 
平均 长 度 Average length (nt) 1 014.54 1047.8 983.83 
中 位 长 度 Median length (nt) 537 574 505 
GC 百分比 GC percentage 42. 7896 42. 1496 41. 9696 
N50 (nt) 1767 1 821 1755 





N50 ; 将 基因 从 长 到 短 排序 ,依次 累加 基因 的 碱 基 数 , 当 累 加 的 碱 基 数 达到 总 碱 基数 50 96 时 最 后 一 个 加 上 的 基因 的 长 度 Length of the last added 


gene when the accumulated base length reached 5096 of the total base length by sorting all genes from the largest to the smallest and accumulating the base 


numbers of the genes in turn. 
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图 1 小 胸 浆 甲 转录 组 基因 与 NCBI 的 NR 数据 库 比 对 结果 
Fig. 1 Alignment of the unigenes from Microdera punctipennis transcriptome against NCBI non-redundant protein sequence database 


A; E fii) fg E-value distribution. 每 个 基因 在 NR 库 中 比 对 得 到 的 最 小 EE 值 序列 为 同 源 序列 (E-value € 1.0E 5) , The sequence of the smallest 
E-value for each unigene is represented as the homolog (E-value «1. 0E-5). B: 同 源 序列 物种 分 布 Species distribution of unigene homologs. 图 中 标 


注 出 每 个 物种 所 占 比 例 。Percentage of species distribution is marked on the chart. 
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2.3 hE FH E SR ZBLTIS] S7 4I B) 28 RES 2T 
在 错误 发 现 率 (FDR) 小 于 0.1 条 件 下 ,利用 
DESeq 分 析 得 到 349 个 上 调 表达 和 70 个 下 调 表达 
基因 ,差异 表达 基因 的 表达 倍数 均 三 3. 6, 即 11og, 
(CD/CK) | =3.6, HP CD 是 低温 处 理 组 的 表达 量 ， 


CK 是 室温 对 照 组 的 表达 量 。 以 DESeq 分 析 的 差异 
A 
DESeq up GFOLDz4 
349 471 


表达 倍数 为 基础 ,将 GFOLD JC PF A3 Pr BO BO EL UE OS 
4,553] 471 个 上 调和 S1 个 下 调 表 达 基 因 。 在 下 游 
分 析 中 以 两 种 方法 得 到 的 差异 性 表达 基因 的 并 集 作 
为 分 析 对 象 。 两 种 方法 共 得 到 514 个 上 调 表达 基因 ，， 
其 中 共同 上 调 表达 306 个 (图 2: A) ; 共 得 到 91 个 下 
调 表达 基因 ,其 中 共同 下 调 表达 60 个 (图 2: B). 


DESeq. down GFOLD <4 


70 81 





图 2 DEseq 和 GFOLD 15:41 Pr^] If E P IGI 26 5e de 3 AE DAL RI 4 E] 
Fig. 2 Venn diagram of the differentially expressed unigenes in Microdera punctipennis transcriptome by DESeq and GFOLD analysis 
A 通过 两 种 方法 共 得 到 的 514 个 上 调 表达 基因 ,其 中 共同 上 调 表 达 306 个 In total, 514 up-regulated unigenes are obtained by using the two tools, 
of which 306 unigenes are common in both analyses; B; 两 种 方法 总 共 得 到 91 个 下 调 表达 基因 ,其 中 共同 下 调 表达 60 个 m total, 91 down- 











regulated unigenes are obtained by using the two tools, of which 60 unigenes are common in both analyses. 


2.4 (RAME JA B] E FR E RA LARRA 
GO 功能 分 类 

鉴于 差异 表达 基因 较 多 ,为 了 更 清晰 地 从 小 胸 
鳌 甲 在 低温 下 启动 或 增强 的 分 子 机 制 方面 认识 其 应 
对 低温 的 分 子 机 理 , 利 用 Blast2go 软件 对 低温 响应 
的 514 个 上 调 表 达 基 因 进 行 功能 富 集 分 析 , 涉 及 生 
物 学 过 程 、 分 子 功能 、 细 胞 组 分 3 大 类 ,35 个 小 类 
(图 3)。 富 集 到 生物 学 过 程 的 基因 占 绝 大 多 数 
(95.72% ) ,其 次 为 分 子 功能 (50. 19% ) 不 同 细胞 
组 分 (31. 71% ) 。 在 生物 学 过 程 方面 ,以 代谢 过 程 
为 主 (26.42% ) ,其 次 为 细胞 过 程 (16. 05% ) 和 单 有 
机 体 过 程 (15. 04% ) , 此 外 ,胁迫 响应 (8. 13% ) Æ 
物 调节 (7. 11% ) 免疫 系统 过 程 (5. 49% ) 定位 
(4.2796 ) .多 有 机 体 过 程 (3.86% ) 也 都 占有 相当 多 
的 比例 。 这 些 结果 表明 ,小 胸 鳖 甲 的 低温 响应 是 生 
物 学 过 程 的 全 方位 响应 ,其 中 以 代谢 过 程 的 改变 为 
主 。 在 分 子 功能 方面 , 以 催化 活性 (46. 1296 ) 和 结 
合 活 性 (39. 53% ) 最 为 显著 ,两 者 共 占 85. 65% ,其 
次 为 结构 分 子 活性 、 分 子 功能 调节 和 转运 活性 , 占 
12.0196 。 在 细胞 组 分 方面 , 富 集 到 细胞 的 最 多 
(31.25% ) , 其 次 为 胞 外 区 域 (23. 3196) 、 细 胞 器 
(18. 496) , 4l HE IR ( 12. 26% )、 大 分 子 复合 物 


(9.8296 ) 等 部 位 (图 3) 。 
2.5 ”低温 胁迫 后 小 胸 整 甲 转 录 组 上 调 表 达 基 因 的 
KEGG 代谢 通路 分 析 

为 了 获得 小 胸 鳖 甲 在 低温 下 启动 或 增强 的 代谢 
通路 ,根据 各 序列 GO 注释 对 应 的 酶 代码 (enzyme 
code, EC) 分 析 小 胸 整 甲 转 录 组 低温 响应 上 调 表 达 
基因 富 集 到 不 同 代谢 通路 的 数目 ,得 到 93 个 基因 富 
集 到 18 个 代谢 通路 (图 4) ,其 中 富 集 前 4 位 的 通路 
是 味 叭 代谢 、 硫 胺 素 代 谢 、 抗 生 素 合 成 、 糖 酵 解 / 糖 异 
生 , 基 因数 目 均 含 10 个 以 上 ; 富 集 前 5 -10 位 的 通 
路 是 氨基 糖 和 核 苷 酸 糖 代谢 、 淀 粉 和 芒 糖 代谢 .果糖 
和 甘露 糖 代谢 、 丁 苷 菌 素 和 新 霉 素 合成 . 半 乳 糖 代 
谢 。 富 集 到 其 他 代谢 途径 的 基因 数目 均 少 于 4 个。 
2.6 小 胸 警 甲 低温 响应 转录 组 中 与 非 生 物 胁 迫 相 

在 胁迫 响应 相关 GO term 中 ,经 Fisher 精确 检 
验 , 共 得 到 8 个 与 非 生物 胁迫 相关 的 显著 上 调 表达 
基因 ,分 别 为 热 激 蛋白 90(HSP90)、 生 长 阻 滞 与 
DNA 损伤 诱导 蛋白 45( Cadd45y) 、 线 粒 体 吡咯 琳 -5- 
羧 酸 脱氧 酶 (P5CDH) 、 超 氧化 物 皮 化 酶 ( superoxide 
dismutase, SOD) 、 钉 联接 蛋白 (calnexin) 和 烯 醇化 酶 
(enolase) 基 因 等 ,上 调 表达 倍数 均 大 于 4.7( 表 2)。 
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代谢 过 程 Metabolic process 

细胞 过 程 Cellular process 

单 有 机 体 过 程 Single-organism process 

应 激 反 应 Response to stimulus 

生物 调节 Biological regulation 

免疫 系统 过 程 Immune system process 

定位 Localization 

多 元 有 机 体 过 程 Multi-organism process 
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图 3 小 胸 警 甲 低温 处 理 转录 组 上 调 表 达 基 因 的 G0 功能 分 类 
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Fig. 3 GO terms alignment of the up-regulated unigenes from the transcriptome of Microdera punctipennis at low temperature 


药物 代谢 -其 他 酶 Drug metabolism - other enzymes 
糖 醇 解 / 糖 异 生 Glycolysis/Gluconeogenesis 

喇叭 代谢 Purine metabolism 
硫 胺 素 代 谢 Thiamine metabolism 
AERE SR E VIG I, Streptomycin biosynthesis 
KARIH Tyrosine metabolism 
抗生素 合成 Biosynthesis of antibiotics 
用 烷 代谢 Methane metabolism 
果糖 和 甘露 糖 代谢 Fructose and mannose metabolism 

光合 生物 碳 固定 Carbon fixation in photosynthetic organisms 

谷 胱 甘 肽 新 陈 代谢 Glutathione metabolism 

泛酸 和 辅酶 A 生物 合成 Pantothenate and CoA biosynthesis 

丁 酰 苷 菌 素 及 新 霉 素 生物 合成 Butirosin and neomycin biosynthesis 
省 粉 和 蔗糖 代谢 Starch and sucrose metabolism 

丙酮 酸 代 谢 Pyruvate metabolism 

半 乳 糖 代谢 Galactose metabolism 

氨基 糖 和 核 苷 酸 糖 代谢 Amino sugar and nucleotide sugar metabolism 
B- 丙 氨 酸 代谢 Beta-alanine metabolism 
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图 4 小 胸 警 甲 低 温 处 理 转录 组 上 调 表达 基因 的 KEGG 38984) (P 10 7^) 


Fig. 4 KEGG pathway analysis of the up-regulated unigenes in transcriptome of Microdera punctipennis at low temperature (P<10™) 
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R2 小 胸 浆 甲 转录 组 中 非 生物 胁迫 响应 GO 类 群 上 调 表 达 基 因 信 息 


Table 2  Unigenes in the abiotic stress response GO term in Microdera punctipennis transcriptome 






































序列 ID E 值 GFOLD fH 序列 描述 

Sequence ID E-value GFOLD value Description 
TR156191c0 gl il 1.5E-25 4. 94623 热 激 和 蛋白 90 Heat shock protein 90 
TR279401|cl gl il 7.9E-114 6.74473 热 激 和 蛋白 90 Heat shock protein 90 
TR279401c0. gl il 0 7.87764 Hsp90 样 蛋白 Hsp90-like protein 
TR25977|cl gl il  6.4E-72 4.7397 HÆRRI- DNA 损伤 蛋白 45 Growth arrest and DNA damage-inducible protein GADD45 gamma 
TR306521c0. gl il 0 5.38352 线粒体 吡咯 琳 -5- 羧 酸 脱氧 酶 Mitochondrial PSCDH 
TR50511c0 gl il 2.1E-66 4.81811 超 氧 化 物 歧化 酶 Superoxide dismutase 
TR293391c0_gl_il 5.4E-112 4.96651 钙 联接 蛋白 Calnexin 
TR79141c0_gl_il 0 5.73647 烯 醇化 酶 Enolase 


3 讨论 


3.1 小 胸 警 甲 4*C 低 温 转录 组 的 质量 分 析 

生活 于 古 尔 班 通 古 特 沙 漠 的 小 胸 浆 甲 适应 于 其 
漠 极 端 温 度 变 化 , 抗 冻 耐寒 性 很 强 ,除了 产生 抗 冻 蛋 
白 渡 过 严冬 之 外 ,体内 还 有 一 系列 生理 生化 变化 , 涉 
及 到 大 量 基因 的 表达 变化 。 由 于 冷 驯化 是 一 种 可 诱 
导 响 应 过 程 ,利用 RNA-sed 技术 不 仅 可 检测 具有 冷 
驯化 效应 的 基因 还 可 检测 其 他 基因 的 表达 变化 ,还 
可 发 现 主 要 控制 冷 驯 化 效应 的 基因 及 响应 一 些 二 级 
变化 的 基因 (Parker et o.，2015 ) 。 本 研究 利用 
Illumina Solexa 测序 平台 对 4% 低温 处 理 3 h 的 小 胸 
浆 甲 成 虫 及 其 25Y 对 照 组 进行 测序 ,得 到 68 165 001 
条 有 效 读 序 , 从 头 组 装 得 到 56 130 条 Trinity 组 装 序 
列 ,利用 CD-HIT 软件 去 除 宛 余 序 列 后 ,最终 得 到 平 
均 长 度 为 984 nt, N50 为 1 755 nt 的 非 宛 余 基 因 
51 712 条 。 同 为 拟 步 甲 科 的 黄粉 虫 触角 转录 组 测序 
得 到 52 216 616 条 有 效 读 序 , 平 均 长 度 451 nt, N50 
为 505 nt 的 基因 35 363 条 (Liu et al., 2015)。 从 获 
得 基因 的 数量 、 平 均 长 度 和 N50 长 度 等 方面 均 显 示 
本 次 小 胸 浆 甲 转录 组 的 拼接 效果 良好 。 研 究 表明 ， 
有 参考 序列 组 装 的 转录 组 敏感 性 高 .了 元 余 度 低 ,并 可 
拼接 出 新 转录 本 (Marchant et al., 2016) 。 基 因 组 已 
测序 的 拟 步 甲 科 的 赤 拟 谷 盗 了， castane 基因 组 大 小 
J 210.265 Mb , 含 18 076 条 和 蛋白 质 ( NCBI)。 与 赤 
拟 谷 盗 基 因 组 信息 比较 ,推测 本 次 拼接 的 小 胸 鳖 甲 
转录 组 可 能 仍 存 在 宛 余 现象 ,但 不 排除 可 变 前 接 的 
可 能 性 。 与 NCBI 的 NR 和 蛋白 库 比 对 可 知 , 小 胸 整 甲 
转录 组 有 29 925 条 基因 与 已 知 蛋 白质 序列 同 源 , 占 
42.13% ,表明 转录 组 中 尚 有 大 量 功 能 未 知 的 新 基 
因 ,这 些 新 基因 可 能 与 昆虫 独特 的 生物 学 特性 和 环 
境 适 应 性 有 关 。 在 已 注释 的 基因 中 , 与 亲缘 关系 最 












































近 的 赤 拟 谷 资 同 源 性 达 72.75% 。 目 前 已 有 23 种 
果 晶 和 21 种 蚊子 等 双 翅 目 昆 虫 的 基因 序列 信息 ,但 
是 鞘 妈 目 昆 虫 的 基因 序列 信息 不 多 ,本 次 高 通 量 测 
序 的 小 胸 鳖 甲 转录 组 数据 将 有 助 于 推进 畏 翅 目 昆 虫 
的 分 子 生物 学 研究 。 
3.2 小 胸 警 甲 低温 转录 组 上 调 表达 基因 的 GO EE 
本 研究 利用 DESeq 分 析 得 到 349 条 上 调和 70 
条 下 调 表达 低温 响应 基因 。 利 用 GFOLD 软件 得 到 
471 条 上 调和 81 条 下 调 表 达 低 温 响应 基因 , 可 见 
GFOLD 软件 的 分 析 结 果 涵 盖 面 更 广泛 。 有 研究 认 
为 GFOLD 软件 在 分 析 非 生物 学 重复 转录 组 时 比 其 
他 和 常用 软件 的 分 析 结 果 更 具 生 物 学 意义 (Feng et 
al., 2012) 。 鉴 于 转录 组 分 析 力 图 获得 更 多 差异 表 
达 基 因 ,我 们 取 两 种 方法 得 到 的 差异 性 表达 基因 的 
并 集 为 分 析 对 象 , 共 得 到 514 个 上 调 、91 个 下 调 表 
达 低 温 响 应 基因 。 本 转录 组 响应 低温 上 调 表达 基因 
的 qRT-PCR 验证 已 在 前 期 的 部 分 工作 中 发 表 ( 李 洁 
HÆ, 2013; Lu et al., 2014a; DEZ RE, 2015; i 
闪闪 等 , 2015 ) 。 对 这 些 上 调 表达 基因 进行 GO 分 
析 , 可 富 集 到 3 大 类 ,35 个 小 类 (图 3)。 其 中 富 集 
到 生物 学 过 程 类 群 占 绝 大 多 数 (95.72% ) ,又 以 代 
谢 过 程 所 占 比例 最 高 ,其 次 为 细胞 内 过 程 类 群 和 单 
有 机 体 过 程 类 群 。 这 与 昆虫 在 低温 下 首先 是 代谢 过 
程 发 生 改 变 是 一 致 的 (Zhang et al., 2011)。 对 寄生 
蜂 Aphidius colemani 蛋白 质 组 学 研究 显示 ,温度 变动 
可 上 调 糖 醇 解 三 羧 酸 循环 `.ATP 合成 相关 的 代谢 通 
路 (Colinet et al., 2007) 。 脂 肪 和 碳水 化 合 物 代谢 与 
冷 暴 露 的 程度 相关 (Zhang et al., 2011), /) NJ E FH 
低温 转录 组 中 胁迫 响应 ` 生 物 调节 和 免疫 系统 过 程 
也 都 占有 重要 比例 。 在 双 翅 目的 果 蝇 中 也 发 现在 低 
温 下 免疫 系统 的 多 个 基因 上 调 表 达 (Zhang et al., 
2011)。 小 胸 整 甲 在 低温 下 免疫 系统 相关 基因 的 上 
调 表达 提示 昆虫 免疫 与 低温 相关 。Parker 等 (2015 ) 
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研究 5% 处 理 6 d HRI D. montana FI D. virilis ,发 
现 分 别 有 177 条 和 458 条 差异 表达 基因 ,GO 富 集 到 
16 个 生物 学 过 程 。 小 胸 鳌 甲 4Y 低温 转录 组 所 获得 
的 上 调 表 达 基 因数 目 和 富 集 的 GO 2S ESAE T AR 
D. montana 和 D. virilis , XX Je T R 8] | Jf] E FH | 3€ 
机 制 更 为 复杂 还 有 待 深 入 研究 。 

小 胸 整 甲 低温 响应 上 调 表达 基因 所 富 集 的 生物 
学 过 程 与 果 晶 D. montana 和 D. virilis 大 体 一 致 ,但 
也 存在 多 方面 不 同 , 表 现在 小 胸 鳖 甲 存 在 显著 的 单 
有 机 体 过 程 ( 占 15. 04%)、 免 疫 系 统 过 程 ( 占 
5.596 ) 和 细胞 杀伤 过 程 , mR D. montana I D. 
virilis 的 低温 转录 组 数据 缺少 单 有 机 体 过 程 和 细胞 
杀伤 过 程 ,免疫 系统 过 程 也 仪 占 0.25% 。 物 种 间 的 
这 种 显著 差异 值得 深入 研究 。 另 一 方面 , 果 蝇 D. 
montana 和 D. virilis 的 低温 转录 组 中 , 发 育 过 程 类 
群 占有 较 高 比例 (19% ) , mi te 7] fa EE FB , x — 2S BE 
仅 占 1.4% o iNNi EFI Aloe D 2E TE E Uh, 而 果 晶 
为 不 耐寒 昆虫 ,它们 在 生物 学 过 程 类 群 方面 的 显著 
差异 表明 昆虫 对 低温 的 基因 响应 存在 物种 差异 性 ， 
即使 在 同属 的 上 述 两 种 果 蝇 之 间 也 有 很 大 不 同 
( Parker et al., 2015) 。 在 小 胸 整 甲 生物 学 过 程 的 分 
子 功 能 方面 ,催化 活性 和 结合 活性 共 占 85.65% ,是 
十 分 显著 的 功能 类 群 , 这 与 其 他 昆虫 在 低温 下 的 分 
子 功能 分 类 结果 类 似 (Zhang et al., 2011; Zhang et 
al., 2012) , 4'C 处 理 的 植物 类 番 阁 前 Solanum 
lycopersicoides 转录 组 有 2 052 个 下 调和 2 409 个 上 
调 表 达 基 因 (Chen et al., 2015) ,其 GO 分 类 与 小 胸 
浆 甲 转录 组 在 一 般 性 生物 学 过 程 、 分 子 功 能 和 细胞 
组 分 的 主要 类 群 基 本 一 致 ,但 是 多 出 了 诸如 细胞 增 




































































疾病 、 药 物 和 化 学 物质 等 数据 库 的 集合 (Kanehisa 
and Goto, 2000) 。 小 胸 鳖 甲 低温 转录 组 上 调 表 达 的 
基因 可 富 集 到 18 个 KEGG 途径 ,其 中 嗓 叭 代谢 、 硫 
胺 素 代 谢 、 糖 酵 解 / 糖 异 生 抗生素 合成 途径 的 富 集 
程度 最 高 。 其 中 , 味 叭 代谢 和 硫 胺 素 代谢 在 其 他 昆 
虫 低温 转录 组 没有 发 现 (Zhang et al., 2011) ; ET 
代谢 与 能 量 代 谢 有 直接 关系 ,昆虫 在 低温 下 首先 是 
停止 或 降低 与 细胞 分 裂 、 组 织 分 化 和 生长 相关 的 主 
要 能 量 消耗 过 程 (Denlinger and Lee, 2010) 。 在 麻 
kaS. bullata 低温 转录 组 中 糖 酵 解 / 糖 异 生 也 具有 显 
著 富 集 的 KEGG 通路 (Teets et al., 2012) ,该 通路 与 
碳水 化 合 物 代谢 直接 相关 , 提示 在 冷 驯 化 过 程 中 ,有 
氧 代谢 转换 为 糖 异 生 / 糖 酵 解 ,进而 启动 合成 碳水 化 
合 物 类 低温 保护 剂 ( Colinet et al., 2007; Zhang et 
al.，2011 )。 而 硫 胺 素 代 谢 也 与 糖 代 谢 相 关 
( Beltramo et al.,2008) ,可 能 也 参与 碳水 化 合 物 类 低 
温 保护 剂 代谢 。 硫 胺 素 代谢 在 小 胸 获 甲 低温 适应 性 
中 的 作用 有 待 深入 研究 。 在 富 集 前 5 — 10 位 的 通路 
中 还 有 氨基 糖 和 核 音 酸 糖 代谢 ` 洗 粉 和 大 糖 代谢 、 果 
糖 和 甘露 糖 代谢 ,反映 了 在 小 胸 警 甲 低温 响应 过 程 
中 , 短 时 间 内 就 可 启动 大 分 子 碳水 化 合 物 的 降解 和 
小 分 子 糖 类 物质 的 合成 ,从 而 维持 细胞 的 渗透 平衡 。 
TEZLFÉ VJ HR S. crassipalpis ( Michaud and Denlinger, 
2007) 和 黑 腹 果 蝇 (Overgaard et al., 2007) 中 也 观察 
到 了 类 似 情况 ,提示 和 葡萄糖 的 产生 可 能 是 冷 驯化 的 
一 般 性 特征 (Teets et al., 2012) 。 此 外 ,小 胸 鳖 甲 中 
氨基 酸 代 谢 丙酮 酸 代谢 也 都 有 较 高 的 富 集 度 。 在 
PRE S. bullata 低温 转录 组 中 发 现 氨基 酸 代谢 通路 
基因 表达 显著 上 调 , 并 且 几 乎 所 有 的 氨基 酸 含量 都 



























































殖 、 细 胞 连接 、 类 核 .电子 载体 活性 、 蛋 白质 结合 转录 
活性 等 类 群 ,这 些 在 小 胸 上 刺 甲 中 所 缺乏 的 GO 类 群 
有 些 是 植物 所 特有 的 。 而 小 胸 整 甲 低温 转录 组 多 出 
了 分 子 功能 调节 融和 行为 两 个 类 群 ,表明 在 动 植物 
之 间 , 虽 然 在 冷 驯化 响应 的 生物 学 过 程 、 分 子 功能 

细胞 组 分 的 类 群 大 致 相同 ,但 也 有 显著 不 同 。 冷 驯 
化 作为 一 种 逆境 适应 性 策略 ,在 生物 中 存在 同 源 性 状 
(Salt et al., 1961; Denlinger and Lee, 2010) 。 由 于 进 
化 速度 不 同 , 同 源 性 在 不 同 水 平 上 可 能 不 一 致 ,直系 
同 源 选 择 或 漂移 可 导致 贡献 同类 生物 过 程 的 基因 被 

















有 增加 ,推测 可 能 与 冷 损伤 修复 有 关 ( Teets et al., 
2012) 。 值 得 注意 的 是 ,小 胸 鳖 甲 冷 响应 上 调 表 达 基 
因 还 富 集 到 了 抗生素 生物 合成 途径 ,这 也 是 与 其 他 昆 
虫 低温 转录 组 的 不 同 之 处 。 利 用 UniProt 数据 库 比 
对 ,结果 显示 ,除了 基因 TR79141c0_gl_il 和 TR16785 
1c0_g2_i4 Zh, Mt KEGG 抗生素 生物 合成 途 
径 的 其 他 高 度 同 源 序列 基因 均 分 布 于 寄生 虫 中 。 昆 
虫 与 微生物 共生 现象 比较 普遍 (Li et al., 2012) ,推测 
荒漠 昆虫 小 胸 鳌 甲 体内 可 能 有 寄生 或 共生 生物 ,它们 
通过 产生 一 些 抗菌 物质 而 保护 自身 及 其 寄主 在 低温 






















































































增 选 或 淘汰 ,使 得 产生 同 源 性 状 的 基因 可 能 不 同 
(Canestro et al., 2007; McCune and Schimenti , 2012) 。 
3.3 1M E FR IG SER B EVE RA EAS] KEGG 
富 集 

KEGG 数据 库 是 一 个 涉及 基因 组 生物 学 途径 、 








下 不 受 其 他 微生物 的 攻击 。 小 胸 获 甲 与 共生 生物 如 
何在 低温 下 的 协同 作用 也 是 一 个 值得 研究 的 领域 。 
3.4 ”小 胸 整 甲 低 温 转 录 组 中 的 非 生 物 胁迫 响应 
基因 

在 胁迫 响应 相关 CO term 中 ,有 8 个 显著 上 调 
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表达 的 非 生物 胁迫 相关 基因 ,包括 3 个 热 激 蛋白 90 
( Hsp90) 基因 ,以 及 生长 阻 滑 和 DNA 损伤 诱导 和 蛋白 
45 (Cadd45y)、 吡 咯 啉 5- 羧 酸 脱氧 酶 (P5CDh) 、 超 
氧化 物 歧 化 酶 (SOD) 、 钉 联接 蛋白 和 烯 醇化 酶 基因 。 
从 功能 分 析 来 看 ,这 些 基因 可 能 从 分 子 伴侣 、 细 胞 周 
期 阻 滞 线粒体 稳定 性 、 抗 氧化 胁迫 等 方面 应 对 低温 
胁迫 。 热 激 和 蛋白 (heat shock protein, HSP ) 是 生物 响 
应 环境 胁迫 时 体内 合成 的 一 类 蛋白质 ,起 分 子 位 倡 
作用 ,防止 其 他 蛋白 质 错 误 折 琶 和 聚集 ,还 参与 蛋白 
质 跨 膜 运输 。 小 胸 鳌 甲 在 4% 处 理 3 h 的 转录 组 中 
有 3 个 Hsp90 基因 显著 上 调 表 达 (GFOLD 值 > 
4.9) ,表明 其 Hsp90 对 冷 驯 化 温度 很 敏感 。 低 温 刺 
激 时 ,不 同 物种 Hsp 的 响应 方式 有 较 大 差异 , 黑 腹 
HE -4% 处 理 30 min 后 ,其 Hsp90 无 差异 性 表 
达 , 而 安第斯 果 量 D. gaucha 的 Hsp90 上 调 表 达 
(Boher et al., 2012) 。 福 寿 螺 Pomacea canaliculata 
在 低温 下 Hsp70 的 表达 没有 发 生变 化 (Zheng et al., 
2012) , [H 1& gj Dl. Placopecten magellanicus 的 低温 处 
理 引起 Hsp70 的 显著 上 调 表达 (Brun et al., 2008) 。 
我 们 在 小 胸 鳌 甲 低 温 响应 转录 组 中 还 发 现 4 条 
Hsp70 转录 本 上 调 表 达 ( GFOLD {Ë > >5.7) ,表明 小 
Ws EE FF A9 Hsp70 和 Hsp90 基因 都 响应 低温 胁迫 。 
IKRE H (calnexin) 也 是 一 种 分 子 伴侣 ,在 内 质 网 
中 与 钙 网 蛋白 (calreticulin ) 一 起 参与 质量 控制 过 
程 , 即 只 允许 正确 折 大 和 组 装 的 蛋白 质 通过 内 质 网 
的 分 泌 通 道 。 小 胸 整 甲 转录 组 中 钙 联 接 蛋 白 基因 呈 
显著 性 上 调 表 达 ( GFOLD ffi =4.966) , X B] t 
甲 在 低温 下 能 够 对 合成 的 蛋白 质 进行 精细 质量 调 
控 , 从 而 保证 新 合成 蛋白 质 的 质量 ( Hammond et al., 
1994)。 生 长 阻 滞 与 DNA 损伤 诱导 和 蛋白 45y 
(Gadd45y) 主要 参与 细胞 周期 阻 滞 、DNA 损伤 修 
复 、 促 进 细 胞 凋 亡 等 功能 (Lucas et al., 2015) 。 其 在 
小 胸 获 甲 低温 转录 组 中 显著 上 调 表达 ,与 低温 下 昆 
虫 的 生长 活动 受 抑制 现象 是 一 致 的 。 昆 虫 在 低温 下 
也 存在 细胞 凋 亡 现象 (Yi et al., 2007) , 7] Ifa Et FH jl 
和 否 存在 细胞 凋 亡 Gadd45y 在 低温 下 的 具体 功能 

竺 深入 人 研究。 此 外 ,昆虫 在 低温 下 ,能 量 代谢 受到 抑 
制 ,虽然 吡咯 啉 5S- 羧 酸 脱氧 酶 P5CDh 的 主要 功能 是 
催化 且 氨 酸 降解 的 第 2 步 以 及 降解 精 氮 酸 为 谷 氮 酸 
( Dendinger and Brill, 1970) ,但 其 在 维持 果 晶 线 粒 
体形 态 及 其 内 环境 稳定 性 方面 的 功能 (He and 
DiMario, 2011) , £z TE fI 3 F7] NA] EE FH PSCDh 可 
能 参与 线粒体 稳定 性 。P5CDh 的 表达 响应 低温 的 
研究 尚未 见报 道 ,该 基因 与 低温 响应 的 关系 及 其 功 
































































































































待 深入 研究。 此 外 ,小 胸 整 甲 低温 转录 组 中 
SOD 基因 的 上 调 表达 ,表明 其 耐寒 机 制 也 包含 抗 氧 
化 胁迫 。 由 于 昆虫 在 低温 下 代谢 活动 降低 ,细胞 内 
会 产生 大 量 活性 氧 分 子 (reactive oxygen species, 
ROS) ,这 些 ROS 分 子 会 对 生物 大 分 子 造 成 氧化 损 
伤 ( Dalton et al., 2009)。 超 氧化 物 歧化 酶 (SOD) 
能 够 催化 超 氧 化 物 自由 基 形 成 过 氧化 氧 ,后 者 再 由 
过 氧化 氢 酶 和 过 氧化 物 酶 催化 产生 水 和 氧气 
(Fridovich, 1977) ,从 而 解除 对 生物 大 分 子 的 过 氧 
化 胁迫 。 小 胸 整 甲 成 虫 低温 响应 转录 组 中 ,SOD 基 
因 显著 上 调 表达 ( GFOLD 值 =4. 818) ,表明 抗 氧化 
是 小 胸 整 甲 低温 冷 驯化 响应 的 机 制 之 一 。 我 们 在 前 
期 研究 发 现 ,小 胸 遇 甲 在 4C 及 -4% 处 理 时 成 虫 体 
内 SOD 的 mRNA 水 平 上 升 ,外 源 表 达 小 胸 整 甲 SOD 
在 低温 条 件 下 可 赋予 大 肠 杆菌 Escherichia coli RE 
酒 酵母 Saccharomyces cerevisiae 细胞 一 定 的 抗 冻 性 
(HR, 2014) 。 果 蝇 在 低温 胁迫 时 SOD 基因 也 显 
著 性 上 调 表 达 并 且 SOD 酶 活性 增高 ( Gao et al., 
2013; Wang et al., 2014) ,这 与 小 胸 效 甲 是 一 致 的 。 
男 外 , 烯 醇化 酶 (enolase) 是 一 种 多 功能 蛋白 ,除了 
作为 糖 酵 解 酶 之 外 ,还 是 纤维 蛋白 溶 酶 原 结 合 蛋 白 ， 
参与 寄生 虫 在 宿主 体内 的 侵 染 和 转移 过 程 。 寄 生 虫 
的 卵 奖 膜 可 产生 烯 醇化 酶 ,在 卵 周围 富 集 纤 维 蛋 白 
溶 酶 原 ,破坏 寄主 胞 外 基质 并 抑制 宿主 免疫 功能 , 协 
助 寄生 虫 逃 避 宿 主 的 免疫 攻击 〈Falabella et al., 
2009; Nguyen et al., 2013) 。 小 胸 鳌 甲 低温 响应 转 
录 组 数据 中 烯 醇化 酶 基因 显著 上 调 表达 ( GFOLD 
值 =5.73647) 。 通 过 GO 富 集 分 析 结 果 可 知 该 序列 
富 集 到 细胞 组 分 大 类 的 卵 绒毛 膜 类 群 。 与 NCBI 和 
SWISS-PROT 数据 库 比 对 结果 显示 , 小 胸 鳖 甲 转录 
组 低温 响应 上 调 表达 的 烯 醇化 酶 可 能 来 自 其 体内 的 
FERo 

本 研究 应 用 RNA-Seq RX ELE Eb a FH 
进行 了 高 通 量 转录 组 测序 分 析 , 共 得 到 68 165 001 有 
效 条 读 序 , 从 头 组 装 获得 平均 长 度 为 984 bp 的 
51 712 条 非 兄 余 基 因 。 其 中 29 925 & ( 24 57. 8796 ) 
与 已 知 蛋 白质 有 同 源 性 。514 条 低温 响应 上 调 表达 
基因 富 集 到 35 个 GO 类 群 ,18 个 KEGG 代谢 通路 。 
其 中 显著 的 GO 类 群 及 KEGG 通路 与 其 他 昆虫 及 植 
物 相 比 都 有 一 定 程度 的 差异 。 在 非 生 物 胁迫 相关 
GO 类 和 群 中 ,有 8 个 基因 显著 上 调 表达 ,包括 3 个 热 
激 蛋 白 90( Hsp90 ) 基因 ,以 及 生长 阻 清和 DNA. 损伤 
诱导 和 蛋白 45( Gadd45y)、 吡 咯 啉 5- 凑 酸 脱 氧 酶 
(P5CDh) JE $E VJ IE f B CSOD ) , 85 HX 2 ÉLRUR 
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醇化 酶 基因 。 表 明 小 胸 整 甲 可 能 从 分 子 伴 倡 .细胞 
周期 阻 灌 线粒体 稳定 性 、 抗 氧化 胁迫 等 方面 应 对 低 
温 胁 迫 。 该 结果 为 描绘 小 胸 警 甲 的 低温 转录 组 图 谱 
及 进一步 探析 小 胸 整 甲 的 耐寒 性 提供 了 基础 ,为 深 
入 研究 范 漠 昆虫 的 耐寒 性 分 子 机 制 及 挖掘 潜在 的 新 
基因 提供 了 丰富 数据 ,对 进一步 研究 该 物种 响应 低 
温 胁 迫 机 制 可 提供 新 的 见解 。 
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